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^ ! ABSTRACT 

Context. The binary fraction among extreme horizontal branch (EHB) stars in Galactic globular clusters (GCs) is an order of magni- 
tude lower than the binary fraction among their counterparts, field hot subdwarfs. This casts serious doubt on their formation channels. 
Aims. In this Letter, I explain the difference between the field and the cluster EHB stars with the binary model of Han et al. (2002, 
2003) for the formation of EHB stars. 

Methods. With the binary population synthesis code of Han et al. (2002, 2003), I follow the evolution of simple stellar populations 
resulting from single star bursts (note that Han et al. 2002, 2003, adopted a constant star formation rate over the last 15 Gyr for the 
i i , production of field EHB stars), and obtain EHB stars at different stellar population ages. 

■ Results. I found that the binary fraction among EHB stars decreases with the stellar population age. The fraction of EHB binaries 
with orbital periods P < 5 d is ~ 2.5% for a stellar population of 10 Gyr from the standard simulation set. 

J«. | Conclusions. The binary model of Han et al. (2002, 2003) is able to explain the lack of EHB binaries in globular clusters. I also 

propose that the precise determination of the physical parameters of close EHB binaries in GCs can lead to the strictest constraint on 

■ common-envelope ejection efficiency. 
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> 1. Introduction 

; m : 

| C*~) Extreme horizontal branch (EHB) stars (also known as hot subd- 
. warfs) are generally considered to be helium-core-burning stars 
' ! with extremely thin hydrogen envelopes (< O.O2M ) (Heber 
t^J- 1986; Saffer et al. 1994). They play an important role in many 
■ ' aspects of astrophysics, e.g., stellar evolution, distance indica- 
' OO , tors, Galactic structure, and the long-standing problem of far- 
, ultraviolet excess in early-type galaxies (Kilkenny et al. 119971 
1 ^ ' Green, Schmidt & Liebert 1986; Han, Podsiadlowski & Lynas- 
; -!h ! Grav l2007l l. 

Various scenarios have been proposed to explain the forma- 
I b ! tion of EHB stars (Mengel, Norris & G ross[l976l Webbink [T984l 
- - 1 Tut ukov & Yungelson [T99Tjl Sweigart [T997l D 'Cruz et al. [l996l 
Lee 119941 Bressan, Chiosi & Fagotto 119941 Yi, Demarque & 
Oemler 1997 ), and most of these scenarios are based on sin- 
gle star evolution (see, however, Mengel, Norris & Gross [1976 ) 
or are meant for single EHB stars (see Webbink 1984; Tutukov 
& Yungelson 1990 for the merger of two helium white dwarfs 
(WDs) to form a single EHB star). Some of the models are also 
quite ad hoc (see Han, Podsiadlowski & Lynas-Gray 2007] for 
a review). Observationally, about two-thirds of field EHB stars 
are in close binaries (Maxted et al. 120011 ). and this presented a 
serious challenge to the theory for the formation of EHB stars. 
Han et al. (2002 2003 ) proposed an a priori binary model (here- 
after HPMM model) for the formation of EHB binaries and sin- 
gle EHB stars. In the model, there are three types of formation 
channels for EHB stars, involving stable Roche lobe overflow 
(RLOF) for EHB binaries with long orbital periods, common- 
envelope (CE) (Paczvnski |1976l ) ejection for EHB binaries with 
short orbital periods, and the merger of helium WDs to form sin- 
gle EHB stars. The model explains the main observational char- 
acteristics of field EHB stars: in particular, their distributions in 



the orbital period-minimum companion mass diagram, and in 
the effective temperature-surface gravity diagram; their distribu- 
tions of orbital period and mass function; their binary fraction 
and the fraction of EHB binaries with WD companions; their 
birth rates; and their space density. The HPMM model is a sig- 
nificant advancement and is widely used in the study of EHB 
stars (O'Tool, Heber & Benjamin 12X3041 ). 

The majority of field EHB stars are in close binaries, and 
HPMM model explained the binary fraction. Therefore, it was a 
great surprise that radial-velocity surveys revealed a remarkable 
lack of close binary systems in EHB stars in globular clusters 
(GCs) (Moni Bidin et al. 120061 Moni Bidin, Catelan & Altmann 
2008). This may imply that the formation channels for cluster 
EHB stars are different from those for field EHB stars. However, 
by checking the channels in the HPMM model and the model 
results, Moni Bidin, Catelan & Altmann (2008) speculated that 
a binary fraction-age relation may be responsible for the lack of 
EHB binaries. Does the relation really exist ? If it exists, is it 
able to explain the lack of EHB binaries in GCs ? 

In this Letter, I carry out detailed Monte Carlo simulations 
for the formation of EHB stars in GCs with the binary population 
synthesis code used for the HPMM model, and obtain the binary 
fraction-age relation, which can explain the lack of EHB binaries 
in GCs. 



2. The model 

In the HPMM model, EHB stars are mainly produced from three 
types of channels: stable RLOF, CE ejection, and the merger 
of helium WD pairs. For a 'primordial' binary system, the pri- 
mary (initially more massive component) evolves and may fill its 
Roche lobe when it has a helium core. Given the mass ratio q of 
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Table 1. Parameters for simulation sets 



Set 


n(q') 








1 


constant 


1.5 


0.75 


0.75 


2 


constant 


1.5 


0.50 


0.50 


3 


constant 


1.5 


1.0 


1.0 


4 


constant 


1.2 


0.75 


0.75 


5 


uncorrelated 


1.5 


0.75 


0.75 



the primary to the secondary less than a critical value, the mass 
transfer is dynamically stable and leads to the formation of a 
wide EHB binary with a main-sequence (MS) companion if the 
helium core is ignited later (1st stable RLOF channel). Given 
the mass ratio over the critical value, a CE is formed, and the 
CE ejection produces a close EHB binary with a MS compan- 
ion if the helium core is ignited later [1st CE ejection channel). 
Some of the 'primordial' binary systems may evolve to WD bi- 
nary systems after the first mass transfer phase. In this case, the 
secondary continues to evolve and may fill its Roche lobe when 
it has a helium core. If the mass transfer is stable, a wide EHB 
binary with a WD companion is resulted given that the helium 
core is ignited later (2nd stable RLOF channel). The mass trans- 
fer is very likely to be dynamically unstable, and this results in 
the formation of a CE, and the CE ejection produces a close EHB 
binary with a WD companion if the helium core is ignited later 
(2nd CE ejection channel). For a WD binary system, if the WD is 
a helium type and the secondary is on the first giant branch (i.e., 
containing a helium core) when it fills its Roche lobe, the resul- 
tant CE contains helium WD pairs. The ejection of the CE leads 
to the formation of close helium WD pairs, and the close helium 
WD pairs may merge due to angular momentum loss of grav- 
itational wave radiation to form single EHB stars (the merger 
channel). 

Binaries with primary's initial mass, Mu ~ 0.95 - 7M and 
initial orbital period Pi ~ 1 - 1000 d, may produce EHB stars 
via the channels described above. For the 1st stable RLOF chan- 
nel, the primary fills its Roche lobe on the Hertzsprung gap or 
the first giant branch, and the mass transfer strips the primary 
of its envelope and leaves a naked helium core. This leads to a 
wide EHB+MS binary, if helium is ignited. For a stellar popula- 
tion of age t < 1 Gyr (Mu ^ 2M ), a binary with P; ~ 1 - 100 d 
can evolve to a wide EHB+MS binary via stable RLOF on the 
Hertzsprung gap or the first giant branch, and a larger EHB 
mass corresponds to a small t (a large Mu). For t > 1 Gyr 
(Mu ^ 2M Q ), only RLOF near the tip of the first giant branch 
(i.e., with a narrow range of P l7 typically Alog(fi/d) ~ 0.5) can 
lead to EHB binaries, as the naked helium core is not ignited if 
the RLOF is not close to the tip of the first giant branch (see 
Table 4 of Han et al. 2002). The corresponding orbital period P\ 
is from ~ 10 d to ~ 1000 d for a stellar population with age t of 
1 Gyr to 15 Gyr. The orbital period and the EHB mass of the re- 
sultant EHB binary would be larger for a large t . For the 1 st CE 
channel to form an EHB binary, the primary of a binary system 
needs to fill its Roche lobe while it has a helium core. However, 
the resultant CE cannot be ejected to form a close EHB binary 
due to a tight envelope if t < 1 Gyr (Mu £ 2M Q ). For t > 1 Gyr, 
the primary of a binary system needs to fill its Roche lobe very 
close to the tip of the first giant branch (i.e., with a typical range 
of initial orbital period Alog(iVd) ~ 0.1) (see Fig. 1 of Han et 
al. |2002| l, otherwise the helium core cannot be ignited after the 
CE ejection. A large t (a small Mu) corresponds to a larger P{ 
(Pi ~ 100 - 1000 d for t ~ 1 - 15 Gyr) and a larger orbital 
period P of the resultant EHB binary. This channel produces a 



close EHB binary with a typical EHB mass of ~ O.46M . No 
EHB star can form from the 2nd stable RLOF channel in the 
HPMM model, as the mass ratio of the secondary (with an ap- 
propriate helium core) to its WD primary is too large and the 
RLOF is not stable. For the second CE channel to form an EHB 
star, the primary of a binary system first needs to experience a 
stable RLOF (with P { ~ 10 - lOOOd) to form a WD binary, 
and the secondary of the WD binary needs to fill its Roche lobe 
when its helium core mass is in an appropriate range (see Fig. 1 
of Han et al.[2002). The ejection of the resultant CE leads to a 
close WD+EHB binary. As the WD can spiral in deeper in the 
envelope during the CE ejection, the WD+EHB can have a much 
shorter orbital period (as short as ~ 0.02 d for a small t) than the 
MS+EHB binary from the first CE ejection channel. For a large 
f, the WD spirals in the envelope of a less massive secondary (the 
envelope is more loosely bound), and the resultant WD+EHB 
binary is wider. For the merger channel, a binary system first 
needs to experience a stable RLOF (with P\ ~ A - 250 d and 
Mu ~ 0.95 - 2 M Q ) to produce a helium WD binary, and the bi- 
nary experiences a CE evolution to form a helium WD pair. The 
pair may coalesce to form a single EHB star due to angular mo- 
mentum loss of gravitational wave radiation. The EHB star pro- 
duced from this channel has a wider mass range (0.4 - O.8M ) 
than from other channels. The mass range is smaller for a small 
t (0.56 - O.64M for t = 2 Gyr) as only very close helium WD 
pairs have time enough to merge, but the range becomes larger 
with a large t (~ 0.4 - O.8M for t = 15 Gyr). 

In the process of CE ejection, the orbital energy released by 
the orbital decay of the embedded binary is used to overcome 
the binding energy of the CE. As is usual, I defined two param- 
eters: the CE ejection efficiency cuce, i.e. the fraction of the re- 
leased orbital energy used to overcome the binding energy; and 
a t h, which defines the fraction of the thermal energy contribut- 
ing to the binding energy of the CE. As in the HPMM model, I 
adopted q crlt -1.5 (for stable RLOF on the first giant branch), 
and cuce = «th = 0.75 as the best choices, and varied them to see 
their effects. 

To obtain the distributions of properties of EHB stars at 
different ages, I have performed detailed Monte Carlo simula- 
tions with the binary population synthesis code developed for 
the HPMM model. In the simulation, I followed the evolution of 
10 million sample binaries according to grids of stellar models 
of solar metallicity and the evolution channels leading to EHB 
stars. I adopted the following input for the simulations (see Han, 
Podsiadlowski & Eggleton |1995l l: 

(1) A single star burst is adopted (note that the star-formation 
rate is taken to be constant over the last 15 Gyr in the HPMM 
model). 

(2) The initial mass function of Miller & Scalo ( 119791 1 is 
adopted. 

(3) For the initial mass ratio distribution n(q') (where q' is the 
ratio of secondary to primary), I adopted (a) a constant mass 
ratio distribution, or alternatively (b) a mass ratio distribution 
where the masses are uncorrelated and drawn independently 
from a Miller & Scalo ( 1979) initial mass function. 

(4) I take the distribution of separations to be constant in log a 
for wide binaries, where a is the orbital separation. The 
adopted distribution implies that ~ 50 % of stellar systems 
are binary systems with orbital periods less than 100 yr. 

Similar to the HPMM model, I have carried out 5 sets of 
Monte Carlo simulations altogether for Population I by varying 
the model parameters over a reasonable range (see Table [TJ. As 
in the HPMM model, the standard simulation set (set 1) has a 
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1st stable RLOF 
1st CE ejection 
2nd stable RLOF 
2nd CE ejection 
merger of HeWD + HeWD 




Fig. 1. The evolution of the fraction f c (in percentage) of EHB 
stars originated from each channel for the standard simulation 
set (set 1), i.e., the 1st stable RLOF channel (wide EHB+MS 
binaries); the 1st CE ejection channel (close EHB+MS binaries); 
the 2nd stable RLOF channel (wide EHB+WD binaries); the 2nd 
CE ejection channel (close EHB+WD binaries); and the merger 
channel (single EHB stars). No selection effect has been applied 
to the lines. Note that there are no EHB binaries from the 2nd 
stable RLOF channel. 




constant initial mass ratio distribution, q CI n = 1.5, and oqe = 
a± = 0.75. 



3. Results 

The simulations give EHB stars at various stellar population 
ages. In order to see the importance of individual evolution chan- 
nels leading to the formation of EHB stars, I plotted in Fig.[T]the 
fractions of EHB stars from different channels (the 1st stable 
RLOF channel for wide EHB+MS binaries; the 1st CE ejection 
channel for close EHB+MS binaries; the 2nd stable RLOF chan- 
nel for wide EHB+WD binaries; the 2nd CE ejection channel 
for close EHB+WD binaries; and the merger channel for single 
EHB stars) at stellar population age t for the standard simulation 
set (set 1). 

Figure [2] shows the evolution of the distribution of orbital 
periods of EHB binaries with stellar population age t for the 
standard simulation set (set 1). The EHB binaries with orbital 
periods P < 5 d can be detected observationally in a GC (Moni 
Bidin, Catelan & Altmann 2008), and I, therefore, showed, in 
Fig. [3j the fractions of the binaries among all the EHB stars, 
including both binaries and singles, at stellar population age t. 
In order to see how the fractions can be affected by model pa- 
rameters, Fig.|3]also displays other simulation sets with various 
model parameters. In the figure, I have applied the GK selection 
effect, which is the selection against EHB stars with compan- 
ions of spectral type G and K (usually MS stars), and is the most 
important selection effect in observations of EHB starfl 



1 We applied the GK selection effect in the following way. If an 
EHB binary has a MS companion and the effective temperature of the 
companion is above 4000 K or the companion is brighter than the EHB 
star, the system is excluded. See the HPMM model for details. 



log (P/d) 

Fig. 2. The evolution of the distributions of orbital periods with 
stellar population age t for the standard simulation set (set 1). 
P is orbital period, and Motai is the total number of EHB stars, 
including both binaries and singles. The stellar population ages 
are denoted and each distribution has been applied an offset of 
-0.1 x (f/Gyr), where t is the age. The distributions are without 
any selection effect. If we apply the GK selection effect, the far- 
right peaks, which are mainly EHB binaries with more massive 
MS companions from the 1st stable RLOF channel, disappear. 
The vertical dotted line denotes P — 5 d. 



o 



■s 



\\y- 




<3*>t 1 


i i \ 




set 2 


i v \ 




. . . ._ set 3 


> ■ 




set 4 


■ \:\ ■ 

!\; \ \. 
1 \'. k 
i V'-.\ 




set 5 


i \" '\ 

; 






\ \ 

i V "• 
\ \ 
V A. 
\ ^ 
\ 







5 10 15 20 



i(Gyr) 

Fig. 3. The evolution of the fraction / (in percentage) of close 
EHB binaries (with orbital periods P < 5 d) among all EHB 
stars, including both binaries and singles. Note that the GK se- 
lection effect is applied. If the selection effect is not considered, 
/ would be smaller due to the fact that the GK selection effect ex- 
cludes wide EHB+MS binaries from the first stable RLOF chan- 
nel. 
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4. Discussion 

Figure [1] shows the fraction of EHB stars originated from each 
channel. We see that the 1st stable RLOF is dominant for t < 
8 Gyr. The resultant wide MS+EHB binaries, however, may not 
be observed due to the GK selection effect as the MS secondaries 
are too massive and bright. The 1st CE ejection channel and the 
2nd CE ejection channel have a contribution peak at t ~ 1.2 Gyr 
and t ~ 3.5 Gyr, respectively. The merger channel starts to con- 
tribute noticeably at t ~ 2 Gyr and dominates at t ^ 8 Gyr, or at 
t 'k. 3.5 Gyr if the GK selection is considered. 

Figure [2] shows the distributions of orbital periods of EHB 
binaries at various ages of a stellar population. For t = 0.5 Gyr, 
the distribution has 3 peaks. The left peak (the small one at 
log(P/d) ~ -1.3) is from the 2nd CE ejection channel, where 
a WD has spiralled in deeply into a tightly bound envelope of 
its companion (with > 2M ) to form a close WD+EHB bi- 
nary. Both the middle peak and the right peak are from the 1 st 
stable RLOF channel with My > 2M , while the middle peak is 
from stable RLOF on the Hertzsprung gap (with a mass ratio of 
qi ~ 1 — 3) and the right peak from stable RLOF on the first giant 
branch (with a mass ratio of q-, ~ 1 - 1.5). For t = 2.5 Gyr, the 
distribution has two peaks. The left peak is from the 1st CE ejec- 
tion channel and the part to the left of the peak is from the 2nd 
CE ejection channel. The right peak is from the 1st stable RLOF 
channel, and the orbital period corresponding to the right peak is 
significantly larger than that at t — 0.5 Gyr, which is due to the 
fact that stellar radius at the tip of the first giant branch is signif- 
icantly larger for My < 2M Q than for My > 2M Q . For t - 5 Gyr 
or larger t, the 2nd CE ejection channel produces WD+EHB bi- 
naries with orbital periods similar to that of MS+EHB binaries 
from the 1st CE ejection channel, as the mass donors (the pri- 
mary for the 1st CE ejection channel, or the secondary for the 
2nd CE ejection channel) are less massive than 2M , and the en- 
velopes are similarly loosely bound near the tip of the first giant 
branch. For a large t, the radius (and the corresponding P\) at 
the tip of the first giant branch is bigger for the donor, and the 
envelope is also more loosely bound, and therefore the left peak 
moves toward a longer orbital period for a large t. 

The EHB binaries with orbital periods P < 5 d are from the 
1st CE ejection channel or the 2nd CE ejection channel. Figure[3] 
shows the evolution of the fraction / of close EHB binaries (with 
P < 5 d) with stellar population age t. We see that / decreases 
sharply with t, which is largely due to the increasing importance 
of the merger channel with t (see Fig. [T). For the standard sim- 
ulation set (set 1), / is ~ 2.5% for t = 10 Gyr. In order to see 
the effects of model parameters, I performed 5 sets of simula- 
tions altogether. Set 2 adopted smaller common envelope ejec- 
tion efficiency and the orbital periods of close EHB binaries are, 
therefore, smaller, leading to bigger / for t > 1 Gyr. Set 3, on 
the other hand, adopted bigger efficiency and, therefore, / is 
smaller. Set 4 adopted smaller q c for stable RLOF on the first 
giant branch, and this makes CE more likely, and leads to more 
EHB binaries from CE ejection channels, which results in bigger 
/. Set 5 adopted an uncorrected initial mass ratio distribution, 
which means a large initial mass ratio qi is more likely and as a 
consequence CE is more likely experienced and, therefore, / is 
larger. Note that the parameters adopted for the standard simula- 
tion set (set 1) produce good matches between observation and 
theory for field EHB stars (see HPMM for details). 

There is a very high fraction of close EHB binaries in the 
field (~ 2/3, see Maxted et al. |2001| l, but t here is a lack of 
close EHB binaries in GCs. Moni Bidin et al. (2006) and Moni 
Bidin, Catelan & Altmann (2008) found that the binary frac- 



tion / (for P < 5 d) is very small in NGC 6752, with a most 
likely value of / = 4% and an upper limit of / = 16% at the 
95% confidence level. They speculated that the sharp contrast 
between field EHB stars and GC EHB stars may be due to a pos- 
sible f-t relatiorQ. They analysed the HPMM model and Fig. 7 
of Han, Podsiadlowski and Lynas-Gray (120071 1 and found that 
such a relation is possible^. This Letter showed that the f-t rela- 
tion clearly exists, in support of the speculation of Moni Bidin, 
Catelan & Altmann (2008). For field EHB stars, their progeni- 
tors have a wide spectrum of masses due to the continuous star 
formation rate during the past, and young (massive) progenitors 
are more likely to produce closer EHB binaries, and therefore 
the binary fraction is much higher (see the HPMM model). 

Gratton et al. (2003] 120031 1 derived the average metallicity 
and the age of NGC 6752 to be [Fe/H] = -1.48 + 0.07 and 
13.8 ± 1.1 Gyr, respectively. The standard simulation set (set 1) 
gives the fraction of close EHB binaries with P < 5 d to be 
/ ~ 1.3% at t — 13.8 Gyr, close to the fraction of 4% observed 
by Moni Bidin, Catelan & Altmann (2008) if we consider that / 
is very sensitive to the model parameters, e.g. CE ejection effi- 
ciency, in each simulation set (see Fig. [3). Note the parameters 
adopted for the standard simulation set give a high binary fre- 
quency among field EHB stars (~ 55% after selection effects 
are taken, see Table 2 of Han et al. [2003), consistent with ob- 
servations of field EHB stars. Note also that our model is for 
solar metallicity. For low metallicity, a star has a more tightly- 
bound envelope. Therefore, the EHB binaries from CE ejections 
have shorter orbital periods, and the helium WD pairs would be 
closer and should be more likely to merge to form single EHB 
stars. However, the production of EHB stars are not expected to 
be affected much by metallicity in the HPMM binary model. 

The CE ejection efficiency strongly affects the distribution of 
close EHB binaries (see Fig . [3j . If we can obtain the parameters, 
i.e., the orbital period, the mass of the EHB star, and the mass of 
the companion, of a close EHB binary in a GC of known age and 
metallicity, we may constrain the CE ejection efficiency. From 
the age and metallicity of a GC, we can infer the mass of the 
progenitor of the EHB star. From the mass of the EHB star and 
the mass of the companion, we would get the location of the CE 
on the evolutionary track of the progenitor, and we would then 
know the binding energy of the CE and the separation before the 
CE. The separation after the CE can be easily calculated from the 
binary parameters observed. Then it would be straight forward 
to derive the CE ejection efficiency. I, therefore, propose here 
that precise determination of the physical parameters of close 
EHB binaries in GCs would possibly give the strictest constraint 
on CE ejection efficiency. However, such an approach seems to 
be too remote. A more realistic observational test of the results 
in this Letter would be searching for EHB binaries with orbtial 
periods of ~ 10 - 20 d in a GC (see Fig. [21). 
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2 Napiwotzski et al. (20041 invoked a possible f-t or /-[Fe/H] rela- 
tion to explain the binary fraction difference between their survey for 
EHB stars and previous surveys for field EHB stars. Their sample con- 
tains more thick disk or halo stars. 

3 Moni Bidin, Catelan & Altmann ( 20081 also argued that dynamical 
effects may not contribute to the formation of EHB stars in GCs, as there 
is a lack of radial gradients among EHB stars in NGC 2808 (Bedin et 
al. |2000t and aj Centauri (D'Cruz et al. |2000l 
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